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Abstract 

Most of the theories involving extra dimensions assume that only the gravitational interaction 
can propagate in them. In such approaches, called brane world models, the effective, 4-dimensional, 
Newton's law is modified at short as well as at large distances. Usually, the deformation of Newton's 
law at large distances is parametrized by a Yukawa potential, which arises mainly from theories 

O '. 

, with compactified extra dimensions. In many other models however, the extra dimensions are 



infinite. These approaches lead to a large distance power-law deformation of the gravitational 
newtonian potential V/v(r), namely V(r) = (1 + kb/r b )V^{r), which is less studied in the literature. 
We investigate here the dynamics of a particle in a gravitational quantum well with such a power- 
law deformation. The effects of the deformation on the energy spectrum are discussed. We also 
compare our modified spectrum to the results obtained with the GRANIT experiment, where the 
effects of the Earth's gravitational field on quantum states of ultra cold neutrons moving above a 
mirror are studied. This comparison leads to upper bounds on b and kf,. 
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I. INTRODUCTION 

Since the pioneering works of Kaluza and Klein jlj], which aimed to unify gravitation 
and electromagnetism in the framework of a 5-dimensional model, theories involving extra 
dimensions have been considerably studied in theoretical physics. In particular, brane world 
,„ode. S , Educed in the ei^ties & are a very act.ve fie.d of research. Although there 
is a great variety of them |3J, |4| , they all rely on the common idea that our 4-dimensional 
universe is a subspace (a brane) of a larger, higher dimensional, universe (the bulk). The 
extra dimensions can either be compactified, or infinite, and the number of extra-dimensions, 
as well as the fields living in the bulk, vary from one model to another. One of their most 
interesting feature is that brane world models are good candidates to solve the so-called 
hierarchy problem, that is the unexplained weakness of the gravitational interaction with 
respect to other forces. 

A first suggestion made by Arkani-Hamed et al jlj] in order to solve the hierarchy problem 

with brane world models is that our universe could be a (4 + n)-dimensional manifold, 

with n compactified extra dimensions 0, Q|. Their basic assumption is that, while all 

the standard model fields - gauge and matter - are confined in the usual 4 dimensions, 

only the graviton can propagate freely in the (4 + n) dimensions, eventually leading to a 

very weak "effective" gravitational interaction in the 4-dimensional universe. Thus, if the 

compactification length scale of the extra dimensions is, say, r c , gravity will only become 

relevant with respect to the other forces at scales lower than r c . Brane world models with 

compact extra dimensions, also known as Large Extra Dimensions (LED) models, and their 

possibly observable consequences, have been considerably studied in the literature . As an 

example, let us consider that the universe possesses n LED, compactified on a torus whose 

n radii are equal to r c . Then, for distances r 3> r c , the newtonian potential 

GmM 

V N {r) = — (1) 



between two pointlike particles is modified to the following form 

V LED (r) = V N (r) (l + «T Ar ), r > r c , (2) 

with a = 2n, A = l/r c . For distances smaller than r c however, gravitation is enhanced since 
the newtonian potential becomes jlj] 

V LE D(r)=V N (r)^j , r<r c . (3) 
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Nowadays, precise constraints on the possible values of a and A in the general Yukawa form 
(J2J) have been obtained from various experiments 9L In the particular case of the LED 
approach, these constraints are upper bounds on rv. JlO. 111 . 12], but in general they can be 
applied to various theoretical models predicting a change in the gravitational law Jr], [h} ]. 

Another approach to solve the hierarchy problem has been proposed by Randall and 
Sundrum in Ref. [l^. Their proposal was also based on brane world models, but they 
allowed for the extra dimension to be infinite. As the LED approaches, brane world models 
with infinite extra dimensions lead to a modification of Newton's law. However, instead of 
the Yukawa form (j2J), these modifications at large distances are generally of the form jj] 



V(r) = V N (r) ( 1 + § 



r > A 



(4) 



where A can be considered as a typical length scale at which the correction due to the 
infinite extra dimensions becomes dominant. For r C A however, V(r) oc r~( n+1 ' as in 
Eq. (j3J), even in this case where the n extra dimensions are infinite. It appears that a power- 
law parametrization such as Eq. (jlj) is poorly studied in the literature compared to the 
Yukawa deformation. Consequently, we think that a discussion of its physical consequences 
is of interest. To this purpose, we propose to investigate the effect of brane world corrections 
such as those of Eq. (jlj) on a particular physical system: The gravitational quantum well 
(GQW), i. e., a bound state of a particle in the gravitational potential. 

A recent experiment, called GRANIT, is devoted to the study of quantum states of 
neutrons in the Earth's gravitational field. Roughly speaking, in this experiment, ultra 
cold neutrons are freely moving (bouncing) in the gravitational field above a mirror. This 
particular setup gives rise to a GQW. As a consequence, the energy spectrum of the neutrons 
in the gravitational field's direction is quantized, and the probability of observing a particle at 
a given height will be maximum at the classical turning points h n = E n /mg, for each energy 



E n . This is precisely what is observed. More details can be found in Refs [16|, jl7j]. This 
experiment gives an opportunity to make a confrontation between observation and various 
theoretical models concerning quantum effects in gravity: Noncommutative geometry 18), 



hat 
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existence of an intrinsic minimal length [19] , spin-gravity effects [20, [2JJ] , . . . Let us note t 
a modification of Newton's law of the Yukawa form (J2J) has already been studied in Ref. 
by using the GRANIT data. Here, we will follow a similar way to derive constraints on the 
corrections (HI). 
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Our paper is organized as follows. In Sec. [TT1 we present the basic properties of the 
gravitational quantum well. Then, we compute the effective potential acting on the bouncing 
particle in Sec. IHIl that is the usual newtonian potential plus a term coming from brane 
world corrections. In Sec. II VI we discuss the effects of such corrections on the GQW spectrum 
and we compute constraints on the correction parameters by comparison with the GRANIT 
experiment. Finally, we draw some conclusions in Sec. |VJ 



II. THE GRAVITATIONAL QUANTUM WELL 

Let us consider the case of a nonrelativistic particle of mass m, subject to a gravitational 
field given by g = —ge z , with g a constant. In order to form a GQW, a mirror is placed 
at z = and actsjis an hardcore interaction. This corresponds to the experimental setup 
described in Refs 



16 



17]]. The corresponding potential is V (x) = V (z) with 

V (z) = +oo for z < 

= mgz for z > 0. (5) 

An infinite potential in z = is a very good description of the mirror, at least for the lowest 
eigenstates. The simplest model describing the vertical motion of the bouncing particle at 
the quantum level is then a Schrodinger equation with a linear potential 

H ^(z)=E ^(z), (6) 

with 

#o = 7T~ + m 9 z - ( 7 ) 
2m 

The mirror, located in z = 0, imposes the boundary condition ^(0) = corresponding to the 
infinite potential barrier. This condition causes the energy to be quantized. The solution of 



such a problem is well-known 
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Problem 40], and reads 

E 0n = ~ (Jj-) a n, (8) 



V> n (*) = N n Ai[6z + a n ], (9) 

where 



a n is the n th zero of the regular Airy function Ai(z) (these zeros can be found for example in 
Ref. {24]), and N n = O 1 ^ 2 /\Ai'(a n )\ is the normalization factor. In analogy with the classical 
situation, a particle of energy E n should bounce at a height given by 

^n = = — 11 

mg 6 

In the GQW, the critical heights are thus quantized. This has been checked by the GRANIT 
experiment [251], in which ultra cold neutrons are moving above a mirror in the Earth's 
gravitational field. The experimentally measured critical heights for the first two states are 

hT P = 12.2/im±1.8 syst ±0.7 stat , (12) 
/i c 2 xp = 21.60m ± 2.2 syst ± 0.7 stat , (13) 

while formula (1111) with m = 939 MeV/c 2 and g = 9.81 m/s 2 gives 

hf = 13.7/im, (14) 
hf = 24.00m. (15) 

The theoretical values are located within the error bars, showing the validity of our sim- 
ple Hamiltonian approach. Actually, as the neutrons are removed from any magnetic field 
in the GRANIT experiment, the main corrections that we could miss with the Schrodinger 
equation ((6]) are spin-gravity corrections, due to the fermionic nature of the neutron. Such 
corrections are of the form S(g x p)/2mc 2 [20j], but they are too small to be observable. As 
a consequence of the good agreement between theory and experiment, the energy shifts due 
to eventual new physical mechanisms must be bounded. A similar argument has already 



been used to derive constraints on theories with a noncommutative geometry 18j, or with 



a minimal length uncertainty 19]. As a possible way to derive constraints on new contri- 
butions, namely power-law modifications of the gravitational potential, we propose here to 
study the variation of 5 = h 2 — h\ when these contributions are added. Theoretically, we 
obtain 5 th = 10.3 0m, while the experimental value is 5 exp = 9.4 0m ±5.4 0m. If we assume 
that the potential in Eq. is no longer mgz but mgz + U(z), with U(z) a perturbation 
term encoding eventual new physics, then 5 will be changed into 5 + Q , with 

n = — [^ 2 \u\^ 2 )-(^\u\^ l )\, (16) 

mg 
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the eigenstates being given by formula 0. The comparison with experimental data requires 
that Q satisfies 

|n| < 4.5 jum. (17) 

From this relation, inequalities ruling the parameters of the perturbation term can be de- 
rived, as we will show in Sec. IIVI Let us point out that working with Q is useful because it 
avoids to deal with eventual arbitrary constants which could appear in U(z). 



III. POWER-LAW CORRECTIONS TO NEWTON'S LAW 



As we already mentioned, all the theoretical approaches involving extra dimensions pre- 
dict a modification of Newton's law. At large distances, the LED framework predicts that 
these deviations should be parametrized by the Yukawa potential 

V(r) = -9H^L (i +ffe ^) , r > r c . (18) 



As the previous deformation has already been studied in Ref. 22j, we will rather consider 
in this work the following form for the modified Newton's law, 

, . GmM ( k b \ , r— , . 

V(r) = + r»yN, (19) 

appearing in the long range corrections due to brane world models with infinite extra di- 
mensions. If Vo(r) = —GMm/r is the usual newtonian potential, we can rewrite Eq. (j4j) as 
V[r) = Vo(r) + Vfe(r). In this picture, A = y/\kb\ is a typical length scale at which the correc- 
tion becomes dominant. Eventually, we will apply our results to the GRANIT experiment. 
As it can be computed from Sec. [Til a typical size of this experiment is around 1/9 ~ 6 /zm. 
As, in every case, the modifications of Newton's law due to extra dimensions are expected 
to become dominant at length scales well below the micrometer (A < 1/5), we assume that 
our analysis is the domain of validity of formula ( |T9l) . 



A. Earth-particle interactions 

Potential (|T9l) describes an interaction between two pointlike particles. However, in our 
case, we are dealing with a pointlike particle bouncing on a plane mirror, at the surface of 
the Earth. Consequently, we have to derive the effective potentials between our particle and 

6 




FIG. 1: Schematic view of the particle, located at a height h above the Earth's surface. 

the Earth on one side, and between the particle and the mirror on the other side. Although 
the mirror's mass is negligible with respect to the Earth's one, its influence should also be 
taken into account because it is much closer to the particle. Firstly, we begin with the 
Earth-particle interactions. Our Planet can be seen in first approximation as a spherical 
body whose density p E is constant, and whose radius is R = 6378 km. If the potential 
between the particle and an infinitesimal Earth's element dM = p E d 3 x' is given by 

then, the total potential acting on the particle is 

r d^T 

V b (x) = -mGp E k b / (21) 
J E \x — x I 

with E the volume the Earth. It can be deduced from Fig. (JI]) that, in spherical coordinates, 

r R r 1 dv 

V b (h) = -2-nmGpEh / r 2 dr / — — , (22) 

Jo J-i[r 2 -2r(h + R)u + (h + R) 2 } {b+1)/2 



where h is the altitude above Earth's surface. For all real values of b ^ 1, 2, 3, we obtain 

(23) 



t/^m 3mgk b 
V b [h) - 



{b-3)R-h (b-l)R + h 



2(b - 1)(6 - 2) (6 - 3)R(h + R) [ h b - 2 (h + 2R) b ~ 2 

where the relations M = (A/3)ttR 3 pe and g = GM/R 2 were used. As a test, we can check 
that we recover the correct Newton's law for fc = 1, 

Tr/IN GmM 

V (h) = -— . (24) 
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For 6 = 1,2, 3, the calculation needs a special treatment but the results are still analytical 

3mgk\ 



V(h) 
V 2 (h) 
V 3 (h) 



4R 
3mgk 2 

2R 
3mgk 3 

4R 



n „ h(h + 2R) 1 fh + 2R 
2R — — in 



In 



h + R 
h + 2R s 



h 



h 



2R 



2R 
h + R 

1 In 



h + 2R 



(25) 
(26) 
(27) 



_h(h + 2R) h + R"'\ h 

Let us note that almost all these expressions are singular in h — 0. 

In our GQW, we can assume that hj R -C 1. By dropping the constants and keeping only 
the dominant terms in h/R, we obtain 



V E (h) 
Vf(h) 

V 2 E (h) 



mgk h 
3mgk 



-h\n[ — 
2R \2R 



h 



3mgk2 
2R 



■ In (6h) 
3mgkf) 



1 



(28) 
(29) 

(30) 
(31) 



2(6- l)(b-2)Rh b - 2 ' 

When it was possible, we replaced \a{h/2R) by Yn(9h) with 9 defined by (JTUl) because 
1/9 is the characteristic length scale of a GQW. The difference in both expressions is just 
an irrelevant constant. 



B. Mirror-particle interactions 

The mirror can be seen as a parallelepiped whose density is pm, that we can parametrize 
as follows: — oo < x,y < oo and —L < z < 0. At the scale of the particle, it is indeed 
justified to approximate the mirror by an infinite plane. In this case, it is easier to consider 
the mirror as a disc with an infinite radius. The equivalent of Eq. fl2"2l is then 

r dr 



V b (h) 



r-0 r R* 
-27rmGpMkb I dz 



(32) 



l [r 2 + {h _ z) 2f+m> 

with R* very large but finite in order to avoid divergent integrals (these divergences will 
disappear later). We find 



3mgk\ 
2R 



L\n{h + L) + Mn 



h + L 
h 



(33) 
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V b>2 (h) 



V 2 (h) : 

3mgkb 



3mg/c 2 ( h + L 
-p In 



2R 
1 



-P- 



1 



W- 2 (h + L) b ~ 2 



(34) 



(35) 



2R r (b-l)(b-2) 

with p = pm/pe- The constant terms involving R* have been removed, since they do not 
have any physical meaning. By only keeping the dominant terms in h/L, we have 



V M (h) « pV E {h) 

v x M {h) ~ 3m9k 



2R 



(36) 
(37) 
(38) 



Let us mention that in the GRANIT experiment, L m 10 cm and p sa 1 22J. One sees 
that the formulae for the mirror have exactly the same analytical form as the ones relative 
to the Earth. This is not surprising since in both cases the size of the neutron is much, 
much smaller than the size of the gravitational sources which appear to the neutron locally 
as an infinite disc. 



IV. THE MODIFIED GQW SPECTRUM 



Influence of the corrections 



We discussed in Sec. (JTXJ) how a new physical mechanism, contained in a perturbation 
U(z), could be included in the GQW. Clearly, in our case, we have 

U b (z) = V b E (z) + V b M (z). (39) 

Roughly speaking, Ub>o(z) oc —mgkb/R. Generally in brane world models, fc& is a positive 
number. In this case, Uf,(z) will bring a negative contribution to the energy, and will thus 
decrease the critical heights. This is coherent with the prediction of almost all modified 
theories of gravitation that the effective gravitational interaction is enhanced, especially at 
short scales. We can remark however that recently, a new model, based in particular on 
Gauss-Bonnet gravity, was proposed, in which the gravitational interaction was decreased 



at short scale 27| . In this special theory, the critical heights would be higher than with the 
usual newtonian potential. At least theoretically, the GQW is able to distinguish between 
theories predicting an enhancement or a decrease of the gravitational interaction because, in 



a classical picture, the particle will be able to bounce higher than predicted with the usual 
Newton's law in the second case, and lower in first case. 

The trivial case where b = does not correspond to a modification of the gravitational 
potential, but rather to a change in the value of the gravitational constant G, which is 
replaced by G(l + ko). It can be encountered for example in theories with supersymmetric 



extra dimensions 



3 if ] . Moreover, the calculations can be done analytically in this case as 



an illustration of our discussion. Thanks to the properties of the Airy function, it can be 
checked that 

<*> = (40) 

and consequently, the level h n will be shifted by an amount 

Ak 

Ah n = -^fh n . (41) 

As we argued previously, if k Q is positive (negative), it will decrease (increase) the criti- 
cal heights. Note that Ah n also increases with n: The corrections logically become more 
important for higher energies. 

Finally, let us remark that, as a level h n will be shifted by Ah n = (Ub) /mg, we have 

Ahn oc h/R. (42) 

This ratio should be very small in any earthly GQW experiment. Consequently, such a 
particular setup does not appear to be very sensitive to the brane world corrections we 
consider here. A confirmation of this point will be given in the next section, where we 
compare our results to the GRANIT experiment. 



B. Constraints on the power-law parameters 



Let us define Ub = —kbmgUb = —A b mgUb- Then, rewriting Eq. (fIBT) . we obtain 



fib = A b Q h 



{HUb\ip2)-(Hu b \A) 



(43) 



By combining the definition (143j) and the experimental inequality (1171) . coming from the 
GRANIT results, we are led to the following constraint on the characteristic length scale of 
the Newton's law modification A, 



|A| < 474.5/ 



/im. 



(44) 
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As the spectrum of the unperturbed GQW is well known, fl b can be computed for every 
value of b in the case of a neutron in the Earth's gravitational field. The upper bounds we 
find on |A| are summed up in Tab. [I] For the calculations we used a ratio between the mirror 
and Earth densities equal to p = 1. In the case b 7^ 1 this means that the mirror gives a 
contribution exactly equal to that of the Earth. For the case b = 1 we chose L — 10 cm. 
Clearly, these bounds are not very good, but this was not a priori known; unfortunately, it 
only appears once the calculations are done. Moreover, we find useful to give them since such 
bounds on power-law deformations are less common in literature than bounds on Yukawa 
deformations. Beyond problems of experimental precision, it is somewhat amazing that 
values of A as large as 1 meter would not affect the GQW spectrum in any observable way, 
since the typical length scale of this system is only 1/6 ~ 6 /im with the parameters of the 
GRANIT experiment. The GQW is thus not very sensitive to the brane world corrections 
that we considered, since these corrections come from effective interactions with the mirror 
and the Earth. As shown in Eq. ( )42l) . these are very small because of the large value of R. If 
the experimental precision is increased of a factor 10-^, the upper bounds on A will roughly 
be improved by a factor 10^ b , thus the improvement will be particularly interesting for the 
smaller values of b. 



b 


1 2 3 4 


|A|(m)< 


8.13 10 4 4.34 0.076 0.009 



TABLE I: Constraints on the power-law parameters of Eq. (|19p . The values p = 1 and L = 10 cm 
were used. 

As an illustration, we can mention some brane world models with infinite extra dimensions 
which cause Newton's law to be modified with Eq. ( fl9l) . We begin with the well-known 
Randall- Sundrum II model [l5], which involves one uncompactified extra dimension. It is 



shown in Ref. 



281 ] that the effective newtonian potential on the brane is 



GmM i 



1 + 4A J? /3vrr) r < A 



V RS = { r v ' ' . (45) 

_GrpM (l+2A 2 /3r 2 ) r>>Afi 

Ar is a typical length scale of the Randall- Sundrum model. Intuitively, one could think 
that it has to be very small, but some recent theoretical arguments suggest that it could be 
about 0.1 mm (^ 10~ 3 eV ), without contradiction with the current observational evidences 

11 



291 ] . The order of magnitude of A# is thus still a matter of discussion. The large r limit 



y. We find then A# < 5.3 m, 



13J]: A R < 0.1 mm. We can 



of potential (|45l) corresponds to /C2 = A 2 = 2A|,/3 respective, 
which is considerably worst than the upper bound of Ref. 
notice that it does not exclude the arguments of Ref. Moreover, even if A# were equal 
to 0.1 mm, it would cause Q to be around 2 fm, which seems clearly unobservable, even 
with an improvement of the experimental precision in GRANIT. Let us point out that the 



Gauss-Bonnet brane world model of Ref. 



271 ] enters in the category 6 = 2. 



Another model which is of interest involves a dilaton field in a 5-dimensional bulk, added 

to the usual gravity [30J]. In this approach, the long-range newtonian potential has the form 

GmM ( A 3 n \ , , . 

V D = r»A D (46) 



Ad depending in particular on the dilaton field 30]. 



Another approach was recently proposed, which consider spherically symmetric solutions 



of Einstein's equations with a bulk cosmological constant and n transverse dimensions 32 , 



331 ] . Gravity is then localised on a four dimensional topological defect. The newtonian 
potential has the form 

GmM I" A n+1 (n)l 
V n = — [l + ^^J, r»A(n). (47) 

Here, the characteristic length scale A(n) depends on the number of extra dimensions. If 
n = 3, we see that this last model predicts a modification in r~ 4 of Newton's law. 

For what concerns theories with b > 5, an important theoretical result has to be pointed 
out. Let us assume, as we did in Sec. IHI] that the characteristic size of the brane world 
corrections, denoted as A, is much more smaller that the gravitational quantum well length 
scale, i.e. 1/9. Then, Eqs ( )28l) and ( !36l) will still be valid at very small h with respect to 
the size of the GRANIT experiment. But, for b > 5, Ub(h) is more divergent at small h 
than h~ 2 . In consequence, if Ub>5 is valid at sufficiently small h, as we assumed here, it will 
forbid the existence of bound states in the GQW, because it is well-known that a potential 
more singular than h~ 2 never admits bound states. Even if A was large enough to avoid 
such a collapse of the bound states, it should bring a very large negative contribution to 
the energy levels. This is obviously in contradiction with what is experimentally observed, 
since the theoretical results obtained with the unperturbed gravitational potential are close 
to the experimental values. We are thus led to the upper bound 

b < 5. (48) 
12 



For example, formula (H7|) predicts that this model should be discarded for n > 3. 

Finally, let us mention that even if we focused on brane world models, our calculations 
can in principle be applied to any theory predicting a new physical mechanism which could 
be added in perturbation as a potential of the form Kb/r h . In particular, our upper bound 
(148p is valid for any real number b. 

V. CONCLUSIONS 

We studied the influence of power-law corrections to Newton's law on the gravitational 
quantum well spectrum. These corrections naturally arise in brane world models with infinite 
extra dimensions, and are generally less studied that the Yukawa ones. A particularity of 
the system we considered is that the interactions are not between pointlike particles, but 
between a pointlike particle and the Earth on one side, and between the particle and a mirror 
on the other side. We have thus computed the effective gravitational potential acting on the 
particle, where the brane world corrections are added in perturbation of the usual newtonian 
potential. Logically, the energy spectrum of our model, and consequently the critical heights, 
are affected by the corrections. In particular, the commonly predicted enhancement of the 
gravitational interaction at short distances causes the critical heights to be decreased: In a 
classical picture, the particle can bounce lower if the gravitation is stronger. 

Our results can be compared to those of the GRANIT experiment, which studies ultra 
cold neutrons in a gravitational quantum well. This comparison allows to obtain upper 
bounds on the typical length scale at which the correction becomes predominant in brane 
world models. These bounds are far from being as precise as those derived in the framework 
of LED models, but this could not be guessed before doing the calculations. This is due 
to the fact that the effective potentials are conversely proportional to the Earth's radius: 
The brane world contribution is thus very small in every case. We can conclude that an 
experiment such as GRANIT is nearly insensitive to such corrections, although its typical 
scale is the micrometer. Newton's law deformations should be better studied with pointlike 
bodies at the micrometer scale if one wants to get more precise bounds on the validity of 
the usual gravitational potential. 

Finally, let us point out that the existence of bound states in the GRANIT experiment is 
sufficient to make the following evident proposition: A theory forbidding bound states should 
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be discarded. This leads us to suggest that a physical brane world model cannot predict 
that at long distance, the correction to Newton's law is more singular that r -5 . Indeed, 
the total effective potential acting on the neutron, i.e. Ub><o(h) should be in this case more 
singular that h~ 2 . This criterion could serve as a test to discard unphysical approaches or 
to restrict the domain of validity of existing models. 
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